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ABSTRACT:  While  living  systems  have  developed  highly 
efficient  ways  to  convert  chemical  energy  (e.g.^  ATP  hydrolysis) 
to  mechanical  motion  (e.g._,  movement  of  muscle)^  it  remains  a 
challenge  to  build  muscle-like  biomimetic  systems  to  generate 
mechanical  force  directly  from  chemical  reactions.  Here  we  show 
that  a  continuous  flow  of  reactant  solution  leads  to  by  far  the 
largest  volume  change  to  date  in  autonomous  active  gels  driven 
by  the  Belousov— Zhabotinsky  reaction.  These  results  demon¬ 
strate  that  microfluidics  offers  a  useful  and  facile  experimental 
approach  to  optimize  the  conditions  (e.g.^  fabrication  methods^  count erions^  flow  rates^  concentrations  of  reagents)  for 
chemomechanical  transduction  in  active  materials.  This  work  thus  provides  much  needed  insights  and  methods  for  the 
development  of  chemomechanically  active  systems  based  on  combining  soft  materials  and  microfluidic  systems. 


■  INTRODUCTION 

This  paper  reports  the  use  of  a  continuous  reactant  flow  in  a 
microfluidic  system  to  achieve  giant  oscillating  volume  changes 
in  spontaneously  active  gels  for  chemomechanical  transduction. 
Conversion  of  the  energy  stored  in  chemical  bonds  into 
mechanical  forces  is  an  essential  process  for  life.  Living  systems 
have  evolved  thermodynamically  efficient  and  environmentally 
benign  ways  of  harnessing  chemical  energy  to  produce  motion^ 
for  example^  using  the  energy  released  by  ATP  hydrolysis  to 
power  the  directed  movement  of  muscle  fibers  or  micro¬ 
tubules.^  On  the  other  hand;  current  man-made  systems  (e.g.; 
internal  combustion  engines);  though  able  to  convert  chemical 
energy  into  mechanical  forces  or  movement;^  have  much  lower 
efficiency  than  biological  systems.  ThuS;  it  is  attractive  to 
construct  biomimetic  systems  that  preserve  the  key  aspects  of 
living  systems  for  generating  mechanical  forces  from  chemical 
reactions.  Inspired  by  a  prominent  feature  of  muscle;  softness; 
researchers  have  extensively  explored  gels  as  a  type  of  synthetic 
material  for  mimicking  muscle.^”^  Among  the  active  soft 
materials;  one  notable  example  is  the  spontaneous  oscillatin| 
gels  driven  by  the  Belousov— Zhabotinsky  (BZ)  reaction^ 
reported  first  by  Yoshida.^  Although  it  successfully  integrates 
chemical  oscillation  with  soft  materials  to  give  beating  that 
resembles  the  behavior  of  heart  muscle  cells;  the  Yoshida  gel 
exhibits  relatively  small  volume  change  ('^17%)  during  each 
cycle.  Subsequent  structural  modification  of  the  polymeric 
networks  in  gels  of  this  type  has  resulted  in  only  modest 
improvement.^'^'^  These  results  imply  that;  in  addition  to 
softness;  it  is  necessary  to  mimic  other  features  of  muscle  in 
order  to  optimize  chemomechanical  transduction. 

In  research  on  biomimetic  gel  materials;  an  often  overlooked 
anatomical  feature  is  that  capillaries  play  an  integral  role  in  vivo 
by  providing  biological  fuels  to  and  removing  metabolic  waste 


from  muscle  cells.  ThuS;  to  explore  gel  materials  in  a  muscle¬ 
like  environment;  we  choose  microfluidics  to  simulate  the 
effect  of  the  capillary  for  chemical  transportation.  Several 
groups  have  explored  gels  in  microfluidic  systems.^^”^^  For 
example;  Beebe  reported  the  use  of  gels  as  a  component  for 
self-regulated  flow  control^  and  adaptive  liquid  microlenses;^ 
Kumacheva^^  and  Seiffert^^  used  microfluidic  systems  to 
generate  capsules  of  biopolymer  hydrogels;  Herr  demonstrated 
the  use  of  gels  for  automated  microfluidic  protein  blotting; 
Wu  utilized  hydrogels  to  create  complex;  static  solution 
gradients  in  microfluidic  channels; and  Ismagilov  employed 
nonlinear  dynamic  reactions  in  microfluidics  to  build  a  minimal 
functional  model  of  hemostasis.  Although  microfluidic 
systems  have  found  application  in  studying  the  chemical 
communication^^  of  BZ  droplets  and  chemical  self-organiza¬ 
tion;^^  the  properties  and  potential  of  self-oscillating  gels  in  a 
microfluidic  system  have  yet  to  be  explored  for  chemo¬ 
mechanical  transduction. 

In  this  work;  we  use  a  microfluidic  channel  (Scheme  l)  to 
explore  chemomechanical  transduction  in  BZ  catalyst-contain¬ 
ing  microgels  made  by  photopolymerization  (Figure  l).  We 
observe  that  the  volume  change  during  each  oscillation  under 
optimal;  continuous  flow  exceeds  500%;  more  than  an  order  of 
magnitude  greater  than  the  highest  reported  value  (25%)^^  for 
volume  changes  of  oscillating  gels  without  flow.  Using  a 
microfluidic  system  to  study  gels  containing  two  new  BZ 
catalysts  (l  and  3;  Figure  2B);  we  confirm  that  continuous  flow 
generally  results  in  larger  volume  changes  than  without  flow. 
Besides  easily  tuning  the  flow  rate;  continuous  flow  allows  us  to 
switch  the  counterion  of  the  BZ  catalyst;  thus  revealing  that  the 
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Scheme  1.  Chemomechanical  Oscillation  of  Active  Gels  with 
and  without  Flow  of  Solution 
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Figure  1.  Setup  for  fabricating  the  microgel  and  the  microfluidic 
channel  for  examining  the  volume  change  of  the  microgel. 


nature  and  the  charge  of  the  electrolyte  significantly  affect  the 
volume  change  of  the  oscillating  gel.  This  work  thus  provides 
new  insights  and  methods  for  the  development  of  chemo¬ 
mechanical  systems  based  on  nonlinear  chemical  dynamics  to 
impart  specific  multiple  functionalities  to  novel  soft  materials. 

■  RESULTS  AND  DISCUSSION 

Figure  1  shows  the  photopolymerization  process  for  fabricating 
microgels.  We  employ  rapid  prototyping^^  to  generate  the 
photomasks^  which  are  intended  to  create  disk-shaped 
microgels  with  diameters  of  100—600  fim.  We  place  the 
solutions  containing  the  catalyst  (l^  2,  or  3)^  the  monomer  (4), 
and;  if  necessary^  the  cross-linker  (S)  (Figure  2A)  into  a 
polycarbonate  mold  (Figures  1  and  SI;  Supporting  Informa¬ 
tion)  with  a  height  of  80  /dm  to  help  define  the  height  of  the 
gel.  After  carrying  out  the  photoinitiated  polymerization^  in  an 
ice— water  bath  for  10—30  min;  we  obtain  donut-shaped 
microgels;  resulting  from  the  reaction  and  diffusion  of  the 
reagents.^^  A  subsequent  3-day  dialysis  in  deionized  (Dl)  water 
ensures  the  removal  of  any  unreacted  reagents  from  the 
microgel. 

Using  the  procedure  in  Figure  1,  we  obtain  three  microgels; 
Gel^  Gelb;  and  Gel^  (Figure  2B).  In  the  three  gels;  the 
monomer;  JV-isopropylacrjlamide  (NIPAAm)  (4);  polymerizes 
to  form  poly(NIPAAm)^^  as  the  polymer  backbone.  In  Gel^ 
and  Gelb;  bis(N;N'-methylene-bis-acrylamide)  (S),  which  serves 
as  the  cross-linker;  and  the  BZ  catalysts  1^^  and  2"^  as  pendants 
copolymerize  into  the  poly(NIPAAm)  to  form  polymer 
network  type  I  (Figure  2C).  Catalyst  3,^  which  bears  multiple 
polymerizable  vinyl  groupS;  itself  acts  as  a  cross-linker;  so  that 
Gel^  only  requires  the  monomer  4  and  the  catalyst  3.  ThuS;  the 
architecture  of  the  polymer  network  in  Gel^.  differs  from  those 
in  Gela  and  Gelb;  and  we  term  it  type  II  in  this  work. 
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Figure  2.  (A)  Molecular  structures  of  the  components  used  for 
fabricating  the  oscillating  gels.  (B)  Compositions  of  the  three 
oscillating  gels.  (C)  Schematic  representations  of  the  polymeric 
networks  in  the  two  types  of  gels.  (D)  Fluorescent  images  of  a  type  I 
microgel  (Gel^)  viewed  from  its  top  and  a  3D  projection  obtained 
from  a  z-scan  (scale  bar  is  100  /im).  (E)  Fluorescent  images  of  a  type 
II  microgel  (Gel^)  viewed  from  its  top  and  a  3D  projection  obtained 
from  a  z-scan  (scale  bar  is  100  /im). 


For  type  I  gels;  the  polymerization  rate  of  NIPAAm  (4)  is 
14%  and  23%  within  Gel^  and  Gel^;  respectively.  The  molar 
ratio  of  catalyst  to  NIPAAm  is  0.5%  and  0.6%  in  Gel^  and  Gel^; 
respectively.  For  the  type  II  gel;  Gel^;  the  polymerization  rate  of 
NIPAAm  is  around  11%;  and  the  molar  ratio  of  catalyst  (3)  to 
NIPAAm  is  around  5%;  almost  10  times  as  high  as  in  the  type  I 
gels.  These  differences  originate  from  the  different  structures  of 
the  catalysts  and  their  different  polymer  networks  in  the  gels. 

Since  the  [Ru(bipy)3]^'^  ^"^-based  BZ  catalysts  (l;^^  2;"^  and 
3^)  fluoresce  upon  irradiation  at  488  nm;  we  use  confocal 
fluorescent  microscopy  to  characterize  the  shapes  and  structural 
properties  of  the  microgels.  As  shown  in  typical  examples 
(Figure  2D;E);  both  types  of  microgels  show  weaker 
fluorescence  in  their  centers  than  at  their  edges.  Together 
with  the  3D  projections;  these  results  suggest  that  both  types  of 
microgels  have  donut-like  shapeS;  that  iS;  they  are  thicker  at  the 
edges  than  in  the  centers.  This  unique  geometry  may  also 
contribute  to  the  large  volume  change  observed;  in  agreement 
with  other  works  demonstrating  geometry-enhanced  volume 
changes.^^  Despite  being  fabricated  under  the  same  conditions; 
the  microgels  exhibit  slight  differences  between  batches. 
Additionally;  inhomogeneity  of  the  light  source  for  polymer¬ 
ization  results  in  a  distribution  of  diameters  and  symmetries 
even  within  the  same  batch  of  microgels  (Figure  S2;  Supporting 
Information). 
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The  length  of  each  dimension  of  the  microgels  is  controlled 
to  be  below  600  the  critical  length  scale  of  the  BZ  wave^  to 
guarantee  the  uniform  oscillation  of  the  microgels.  To  minimize 
the  influence  of  size  variation  of  the  microgels  on  their 
oscillating  volume  change^  microgels  with  the  same  size  were 
picked  for  further  tests. 

Fluorescent  images  of  thin  slices  of  the  microgels  (Figure  S3; 
Supporting  Information)  reveal  different  morphologies  of  the 
two  types  of  gels.  While  the  distribution  of  [Ru(bipy)3]^^ 
appears  to  be  homogeneous  in  type  I  microgels  (Gel^  and 
Gelb);  type  II  microgels  (Gel^,)  contain  heterogeneously 
distributed  [Ru(bipy)3]^^  centers^  which  agrees  with  our 
previous  observation.^ 

To  test  the  chemomechanical  oscillation  behavior  of  the 
microgels  in  the  presence  of  floW;  we  use  a  sandwich  type  of 
flow  cell;  constructed  by  gluing  two  pieces  of  polycarbonate 
film  onto  a  250  fim  high  spacer  (Figure  l).  A  syringe  pump 
connected  to  the  inlet  of  the  flow  cell  controls  the  flow  rate  of 
BZ  solution  through  the  cell.  This  simple  setup  (Figure  S2; 
Supporting  Information)  fits  under  an  optical  microscope  to 
allow  an  adapted  CCD  camera  to  record  the  oscillation  of  the 
microgel.  The  flow  cell  is  placed  on  a  temperature  control  stage 
(Figure  S5;  Supporting  Information)  to  maintain  constant 
temperature  during  the  oscillation  of  the  microgel.  Using  this 
microfluidic  system;  we  first  tested  Gel^  under  various  BZ 
reaction  conditions  (Figure  S6;  Supporting  Information)  (e.g.; 
[NaBr03]  =  0.08-0.2  M;  [CH2(C00H)2]  =  0.06-1.2  M;  and 
[HNO3]  =  O.I— 0.8  M);  at  different  temperatures  (Figures  3A 
and  S7;  Supporting  Information).  With  the  same  flow  of  BZ 
reagentS;  Celt’s  oscillating  volume  change  decreases  gradually 
when  the  temperature  rises  and  is  smallest  at  30  °C.  Gel^ 
exhibits  its  maximum  volume  change  during  oscillation  with 
[NaBr03]  =  0.2  M;  [CH2(C00H)2]  =  0.4  M;  and  [HNO3]  = 
0.4  M;  at  10  °C;  different  from  its  optimal  temperature  for 
oscillation  in  static  BZ  solution;  18  °C. 

To  determine  the  optimal  temperature  for  the  other 
oscillating  gels;  Gel^  and  Gel^  were  tested  under  the  same 
continuous  flow  at  different  temperatures.  As  shown  in  Figure 
3B;  Gelb  achieves  its  largest  volume  change  at  22  °C.  Both 
lowering  and  raising  the  temperature  decrease  the  amplitude  of 
oscillatory  volume  change  dramatically.  Like  Gel^;  Gel^  has  a 
type  I  polymer  network  (Figure  2C);  but  the  linker  between  the 
catalyst  and  the  polymer  backbone  is  shorter  than  in  Gel^.^^  In 
Figure  3C;  the  type  II  microgel;  Gel^  also  has  a  maximal 
volume  change  at  10  °C.  However;  the  volume  change  of  Gel^ 
decreases  more  rapidly  with  temperature  rising  than  does  that 
of  Gel^.  Moreover;  there  is  almost  no  mechanical  oscillation 
when  the  temperature  is  higher  than  22  °C.  Similar  to  our 
previous  results  on  type  II  gels;  high  temperature  produces  a 
dense  and  rigid  network  of  Gel^  that  inhibits  the  oscillation. 
Because  the  ratio  of  monomers  to  additives  in  the  polymer  also 
affects  the  temperature  response  of  the  gels;  it  is  expected  that 
the  three  microgels  would  have  different  profiles  of  volume 
changes. 

Comparing  both  types  of  gels'  maximal  volume  change  under 
continuous  flow  of  BZ  solution  to  their  performance  in  static 
solution;  continuous  flow  generally  results  in  larger  volume 
changes.  As  summarized  in  Figure  4A;  under  continuous  floW; 
the  average  diameter  change  of  Gel^  reaches  38%  between  the 
oxidized  state  and  reduced  state;  which  corresponds  to  a  163% 
volume  change.  In  a  bigger  container  filled  with  reaction 
solution  without  floW;  Gel^  exhibits  only  a  4.6%  volume  change 
during  oscillation.  The  continuous  flow  results  in  more  than  an 
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Figure  3.  Relative  volume  changes  of  oscillating  (A)  Gel^;  (B)  Gek; 
and  (C)  Gek  under  continuous  flow  of  BZ  solution  (lO 
[NaBr03]  =  0.2  M,  [CH2(C00H)2]  =  0.4  M,  and  [HNO3]  =  0.4  M) 
at  different  temperature. 


order  of  magnitude  increase  of  the  volume  change.  Under 
continuous  floW;  the  average  diameter  change  of  Gel^  between 
its  oxidized  and  reduced  states  reaches  100%;  corresponding  to 
more  than  a  500%  volume  change.  Without  the  flow  (i.e.;  in 
static  BZ  solution);  Gel^  exhibits  a  25%  volume  change  (the 
best  case).  Unlike  Gel^;  Gel^  has  a  type  II  polymeric  network 
(Figure  2C)  and  shrinks  in  the  oxidized  state.^  Under 
continuous  floW;  the  average  diameter  change  of  Gel^  is  16%; 
which  corresponds  to  a  56%  volume  change.  However;  in 
stationary  BZ  solution;  Gel^,  only  exhibits  a  5%  relative  diameter 
change;  which  corresponds  to  a  16%  volume  change. 

Besides  inducing  large  oscillating  volume  changes;  continu¬ 
ous  flow  slows  the  rate  of  oscillation  of  the  gels.  The  oscillation 
periods  of  Gel^  Gel^;  and  Gel^  are  around  2000;  3600;  and  1500 
s  under  optimal  continuous  floW;  but  they  are  around  800;  250; 
and  400  s,  respectively;  under  static  conditions.  These  results 
further  confirm  that  the  flow  decreases  the  rate  of  accumulation 
of  the  intermediates  of  the  BZ  reaction.  By  examining  the 
correlation  between  oscillating  volume  change  and  oscillation 
period  (Figure  5);  we  also  found  that  larger  volume  change 
corresponds  to  a  longer  oscillation  period  and  relatively  larger 
standard  deviation. 

Considering  the  kinetic  nature  of  the  gel  networks;  it  is  not 
surprising  that  we  observe  batch  differences  and  relatively  large 
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Figure  4.  (A)  Relative  volume  changes  of  oscillating  gels  at  optimal 
rate  of  continuous  flow  and  in  static  solution.  (B)  Volume  change 
ratios  for  gels  between  optimal  continuous  flow  and  static  solution 
largest  volume  change  in  previously  reported  static  results 
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Figure  5.  Temperature-dependent  correlation  between  the  oscillating 
volume  change  and  mechanical  oscillation  period  of  Gel^  Gel^;  and 
Gel,. 


standard  derivations.  Since  the  larger  volume  change  and  longer 
oscillation  period  correlate  with  longer  relaxation  times^  the  gel 
networks  also  evolve  with  the  chemical  oscillation.  This 
intriguing  phenomenon  becomes  even  more  dramatic  when 
the  gel  oscillates  under  continuous  flow^  further  indicating  the 
unique  nature  of  spatially  distributed  chemical  oscillation  under 
flow. 

By  comparing  the  shapes  and  volumes  of  the  microgels  in  the 
reduced  state  at  the  beginning  of  the  BZ  oscillation  and  after 
several  hours  of  oscillation  under  continuous  flow  (Figure  SIO^ 
Supporting  Information)^  we  found  that  the  volumes  of  all 
three  kinds  of  microgels  increase  considerably.  This  observation 
indicates  that  the  networks  of  the  microgels  become  less  rigid 
during  the  BZ  reaction  due  to  the  dynamic  nature  of  the  gels. 


Because  the  microgels  with  larger  oscillation  volume  change 
have  longer  oscillation  periods^  more  extensive  morphological 
changes  occur  in  these  gels.  These  morphological  dynamics^ 
resulting  from  larger  changes  of  their  polymer  networks  in  the 
gel;  thus  contribute  to  the  larger  standard  deviation. 

The  continuous  flow  also  reveals  that  the  geometry  of  the 
microfluidic  channel  influences  the  effect  of  flow  rate  on  the 
volume  change  (Figure  6).  As  shown  in  Figure  6 A,  when  the 


2.5  pL/min  5  pL/min  10  pL/min  50  pL/min 
Flow  Speed 


B  % 


2.5  pL/min  5  pL/min  10  pL/min 

Flow  Speed 


Figure  6.  Relative  diameter  changes  of  Gek  in  a  microfluidic  channel 
with  wid  =  40  (a)  and  a  channel  with  w/d  =  10  (B)  at  different 
continuous  flow  speeds  (w,  width  of  microfluidic  channel;  d,  diameter 
of  microgel;  per  oscillation). 


width  of  the  microfludic  channel  is  much  larger  than  the 
diameter  of  the  microgel  (e.g.;  wjd  =  40),  Gel^  shows  no 
mechanical  oscillations  but  only  chemical  oscillations  at  flow 
speeds  of  2.5  and  5  ^L/min.  The  biggest  volume  change  of 
mechanical  oscillation  happens  at  a  flow  speed  of  10  ^L/min. 
Above  10  fiL/min,  changing  the  flow  speed;  even  by  a  factor  of 
5;  has  little  effect  on  the  volume  change  of  mechanical 
oscillation.  However;  when  the  vndth  of  the  microfludic 
channel  is  relatively  smaller  (e.g.;  w/d  =  lO);  mechanical 
oscillation  occurs  at  a  flow  speed  of  2.5  fiL/ min;  and  the  biggest 
volume  change  is  achieved  at  5  ywL/min.  Compared  to  w/d  = 
40;  when  w/d  =  10;  doubling  the  flow  rate  significantly  affects 
the  magnitude  of  the  volume  change  (Figure  6B).  ThuS;  higher 
flow  speed  is  not  only  necessary  but  also  ideal  for  mechanical 
oscillation  when  the  vndth  of  the  microfluidic  channel  is  much 
larger  than  the  diameter  of  the  microgel.  These  experiments 
demonstrate  that  it  is  the  fluid  velocity  in  the  channel  that  is 
important  and  not  the  flow  rate. 

Continuous  flow  contributes  to  two  processes  that  have 
opposite  effects  on  the  oscillation  of  the  gel.  First;  flow  delivers 
the  reagents  needed  in  the  reaction  [Br03“;  CH2(C00H)2; 
and  HNO3].  If  the  gel  is  operating  in  a  diffusion-dominated 
regime;  then  increasing  the  flow  rate  will  increase  the 
concentration  of  reagents  at  the  gel— solution  interface;  thereby 
increasing  the  rate  of  reaction.  Eventually;  the  rates  will  become 
reaction-limited;  and  increasing  the  flow  beyond  this  point  will 
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not  increase  the  chemical  concentrations.  However^  flow  also 
has  an  antagonistic  effect  on  the  oscillation  through  removal  of 
essential  products  of  the  BZ  reaction^  such  as  Br2;  that  are 
produced  inside  the  gel  and  diffuse  into  the  fluid  stream  and  are 
subsequently  swept  away.  Therefore^  the  reaction  will  cease 
when  the  flow  rate  is  too  high.  With  too  little  flow,,  the  gel  is 
starved;  and  with  too  much  floW;  essential  chemicals  produced 
by  the  gel  are  removed.  ThuS;  it  is  reasonable  that  there  is  an 
intermediate  flow  rate  for  which  there  will  be  a  maximal  volume 
change. 

The  continuous  flow  provided  by  the  microfluidic  setup  also 
allows  quick  and  convenient  screening  or  testing  of  other 
parameters  of  oscillating  gels  for  chemomechanical  conversion. 
For  example;  to  probe  the  role  of  the  counterion  of 
[Ru(bipy)3]^'^;  we  flow  acids  with  anions  of  different  valence; 
in  the  sequences  HNO3— H2SO4  and  H2SO4— HNO3.  That  iS; 
after  the  microgel  (Gel^)  exhibits  relatively  stable  oscillation  in 
HNO3  (0.4  M);  we  change  the  HNO3  to  H2SO4  (0.2  M)  in  the 
flow  stream  of  the  BZ  solution.  Similarly;  after  the  oscillation  of 
the  microgel  becomes  stable  in  H2SO4  (0.2  M);  we  change  the 
acid  back  to  HNO3  (0.4  M).  The  relatively  low  concentrations 
of  both  acids  ensure  the  constant  acidity  of  the  BZ  solution.  As 
shown  in  Figure  lA,  replacement  of  nitric  acid  by  sulfuric  acid 


Figure  7.  Chemomechanical  oscillation  profile  during  transition  of 
Gel^  induced  by  counterion  exchange  (A)  from  nitric  acid  (HNO3)  to 
sulfuric  acid  (H2SO4)  and  (B)  from  H2SO4  to  HNO3. 


biological  phenomena; ranging  from  the  collective  behavior 
exhibited  by  single  cell  organisms  to  neuronal  circuits.  Inspired 
by  thiS;  we  combined  microfabrication  and  microfludics  to 
make  possible  the  study  of  the  oscillatory  behavior  of  an  array 
of  microgels  under  continuous  flow.  As  shown  in  Figure  8A;  we 
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results  in  a  decrease  of  the  volume  change  (e.g.;  from  60%  to 
20%  relative  diameter  change);  accompanied  by  more  rapid 
oscillation  (e.g.;  the  period  decreases  from  2000  to  1300  s, 
Figure  S12A;  Supporting  Information).  This  process  is 
reversible;  though  not  symmetric;  as  switching  sulfuric  acid  to 
nitric  acid  increases  the  volume  change  (e.g.;  from  5%  to  60% 
relative  diameter  change;  Figure  7B)  and  slows  the  oscillation 
(e.g.;  the  period  increases  from  190  to  3500  s,  Figure  S12B; 
Supporting  Information). 

The  smaller  volume  change  caused  by  the  use  of  sulfate 
probably  arises  from  the  buffering  effect  of  the  equilibrium 
between  HS04“  and  S04^“;  which  makes  the  polymer  network 
less  sensitive  to  changes  in  the  charge  state  of  the  catalyst  (e.g.; 
between  [Ru(bipy)3]^'^  and  [Ru(bipy)3]^'^). 

Interestingly;  it  takes  about  4000  s  to  stabilize  the  oscillation 
of  the  gel  when  nitric  acid  is  replaced  by  sulfuric  acid  but  only 
450  s  for  the  reverse  transition  from  sulfuric  acid  to  nitric  acid. 
This  result  indicates  that  the  in-gel  ion  exchange  from  N03“  to 
S04^“  is  much  slower  than  from  S04^“  to  N03“;  which  is  the 
first  observation  to  date  on  the  dynamic  ion  exchange  that 
occurs  in  self-oscillating  gels.  These  results  confirm  that  the 
character  of  the  electrolyte  also  impacts  chemomechanical 
conversion  in  oscillating  gels. 

Coupled  nonlinear  chemical  oscillators  without  flow  have 
been  studied  as  experimental  models  for  understanding 


Figure  8.  (A)  Optical  images  of  an  array  of  Geft  at  different  oscillation 
states  under  continuous  flow  (lO  //L/min).  (B)  Relative  diameter 
changes  of  oscillating  microgels  at  different  positions  in  the  array.  (C) 
Oscillation  profiles  of  the  microgels  (at  positions  I  — IV)  in  the  array. 


made  an  array  of  microgels.  By  analyzing  the  oscillations  of 
microgels  at  positions  I— IV;  we  found  that  microgels  at 
different  positions  in  an  array  oscillated  with  similar  relative 
volume  change  (Figure  8B).  A  plot  of  the  time-dependent 
volume  change  (Figure  8C)  of  the  microgels  indicates  that  the 
microgels  swell  and  deswell  nearly  simultaneously.  The  origin 
of  coupling  between  oscillators  without  flow  is  diffusive  flux  of 
the  participating  components.  It  is  known  that  the  geometric 
arrangement;  size;  and  shape  of  the  oscillators  and  the  distance 
between  oscillators  all  affect  the  coupling  strength  and  the 
resulting  collective  behaviors  of  the  oscillators.^  Applying  a 
continuous  flow  of  reagents  to  an  oscillating  microgel  array 
significantly  modifies  the  diffusive  coupling  of  the  participating 
species.  Depending  on  the  direction  and  speed  of  the  floW;  the 
coupling  strength  between  oscillating  microgels  can  be 
manipulated.  As  the  first  example  of  coupled  oscillating 
microgel  oscillators;  this  system  provides  a  promising 
experimental  model  as  a  tunable  coupling  system  that  deserves 
further  detailed  study. 
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■  CONCLUSION 

This  work;  as  the  first  example  of  using  microfluidics  to  study 
chemomechanical  transduction  based  on  the  BZ  reaction^  sets  a 
new  benchmark  for  the  volume  change  of  self-oscillating  gels. 
Compared  with  classical  thermal  polymerization  methods  for 
the  fabrication  of  BZ  gels;  photopolymerization  has  several 
advantages  (e.g.;  speed;  shape  control;  and  size  control);  which 
should  be  useful  for  further  exploration;  such  as  the  fabrication 
of  hybrid  oscillating  gels  consisting  of  different  types  of 
polymeric  networks.  The  autonomous  large  volume  change 
obtained  in  this  work  also  differs  fundamentally  from  the  shape 
change  of  gels  resulting  from  the  alternating  flow  of  different 
solutions  (e.g.;  buffers  with  different  pHs^^).  It  is  surprising  that 
the  parameters  of  the  continuous  flow  corresponding  to  the 
microenvironment  of  the  oscillating  microgel  affect  so  strongly 
their  chemomechanical  transduction. 

Although  it  is  known  that  the  volume  change  can  be  related 
to  the  size  of  the  BZ  gel  since  traveling  chemical  waves  give  way 
to  uniform  oscillations  in  smaller  systems;^^  the  large  difference 
in  volume  change  observed  in  this  case  is  unlikely  to  be  caused 
by  size  differences  in  the  microgels;  since  the  microgels  used 
have  roughly  the  same  size.  The  fact  that  the  variation  of  the 
flow  speed  has  a  dramatic  influence  on  the  volume  change 
suggests  that  the  difference  in  volume  change  results  primarily 
from  the  floW;  which  alters  the  reaction— diffusion  profiles  of 
the  reactantS;  intermediates;  and  products  of  the  BZ  reaction. 

The  dependence  of  the  volume  changes  of  the  gel  on  the  rate 
of  the  flow  of  the  BZ  reagents  implies  that  gradients  of  the 
reagentS;  intermediates;  and  products  inside/outside  of  the  BZ 
gels  play  an  important  role.  Measurement  of  these  gradients  is  a 
challenging  experimental  problem;  which  lies  beyond  the  scope 
of  the  present  work. 

The  combination  of  microfabrication  and  microfluidic 
techniques  provides  new  and  useful  experimental  tools  for 
the  development  of  theoretical  frameworks  for  understanding 
chemomechanical  materials  based  on  a  wide  range  of  catalysts 
and  polymeric  networks.^'^^ 

■  ASSOCIATED  CONTENT 
O  Supporting  Information 

Synthetic  procedures;  characterization  methods;  Figures  SI  — 
SI 2;  and  two  videos  showing  the  volume  changes  of  Gel^  and 
an  array  of  Gel^.  This  material  is  available  free  of  charge  via  the 
Internet  at  http://pubs.acs.org. 

■  AUTHOR  INFORMATION 
Corresponding  Author 

bxu(^brandeis.edu;  epstein(Sbrandeis.edu;  fraden(^brandeis. 
edu 

Author  Contributions 

^These  authors  contributed  equally. 

Notes 

The  authors  declare  no  competing  financial  interest. 

■  ACKNOWLEDGMENTS 

This  work  was  partially  supported  by  grants  from  the  Army 
Research  Office  (ARO  56735-MS);  the  National  Science 
Foundation  (CHE-1012428);  the  Pioneer  Research  Center 
Program  of  the  National  Research  Foundation  of  Korea  (2012- 
0001255);  and  by  an  NSF  MRSEC  grant  (DMR-0820492).  We 
thank  Prof  Sukho  Park  for  assisting  the  fabrication  of  the  gels 
used  in  Figure  8. 


■  REFERENCES 

(1)  Lodish;  H.;  Berk;  A.;  Kaiser;  C.  A;  Krieger;  M.;  Bretscher;  A.; 
Ploegh;  H.;  AmoU;  A;  Scott;  M.  P.  Molecular  Cellular  Biology,  7th  ed.; 
W.  H.  Freeman:  New  York;  2012. 

(2)  DuaU;  W.;  LiU;  R;  Sen;  A  /.  Am.  Chem.  Soc.  2013;  135;  1280. 
Sengupta;  S.;  Dey  K.  K.;  Muddana;  H.  S.;  Tabouillot;  T.;  Ibek;  M.  E.; 
Butler;  P.  J.;  Sen;  A  /.  Am.  Chem.  Soc.  2013,  13S,  1406.  Paxton;  W.  F.; 
Kistler;  K  C.;  Olmeda;  C.  C.;  Sen;  A;  St.  AngelO;  S.  K.;  CaO;  Y.; 
Mallouk;  T.  E.;  Lammert;  P.  E.;  Crespi;  V.  H.  /.  Am.  Chem.  Soc.  2004, 
126,  13424.  Wang;  Y.;  Fei;  S.-t.;  Byun;  Y.-M.;  Lammert;  P.  E.;  Crespi; 
V.  H.;  Sen;  A;  Mallouk;  T.  E.  /.  Am.  Chem.  Soc.  2009;  131,  9926. 

(3)  Osada;  Y.;  Okuzaki;  H.;  Hori;  H.  Nature  1992,  355;  242.  OterO; 
T.  F.;  Sansinena;  J.  M.  Adv.  Mater.  1998;  10,  491.  Yu;  Y.  L.;  Ikeda;  T. 
Angew.  Chem.-Int.  Ed.  2006;  45,  5416.  Shinohara;  S.;  Seki;  T.;  Sakai;  T.; 
Yoshida;  R;  Takeoka;  Y.  Angew.  Chem.-Int.  Ed.  2008;  47,  9039.  ZhaO; 
Y.  L.;  AprahamiaU;  L;  Trabolsi;  A;  Erina;  N.;  Stoddart;  J.  F.  /.  Am. 
Chem.  Soc.  2008;  130,  6348.  Liang;  S.  M.;  Yu;  Qi  M.;  Yin;  H.  Y.;  Wu;  Z. 
L.;  Kurokawa;  T.;  Gong;  J.  P.  Chem.  Commun.  2009,  7518.  ZhoU;  Y.  X.; 
Sharma;  N.;  Deshmukh;  P.;  Lakhman;  R  K.;  Jain;  M.;  Kasi;  R.  M.  /. 
Am.  Chem.  Soc.  2012;  134,  1630.  Dhanarajan;  A.  P.;  Misra;  G.  P.; 
Siegel;  R  A.  /.  Phys.  Chem.  A  2002,  106,  8835.  Yoshida;  R.  Adv.  Mater. 

2010,  22,  3463.  McDowell;  J.  J.;  Zacharia;  N.  S.;  PuzzO;  D.;  Manners; 
L;  Ozin;  G.  A.  /.  Am.  Chem.  Soc.  2010;  132;  3236.  Zhang;  Y.;  Li;  N.; 
Delgado;  J.;  GaO;  Y.;  Kuang;  Y.;  Fraden;  S.;  Epstein;  1.  R;  Xu;  B. 
Langmuir  2012,  28,  3063.  Suzuki;  D.;  Taniguchi;  H.;  Yoshida;  R.  /.  Am. 
Chem.  Soc.  2009,  131,  12058. 

(4)  Yoshida;  R;  Takahashi;  T.;  Yamaguchi;  T.;  IchijO;  H.  /.  Am.  Chem. 
Soc.  1996;  118,  5134. 

(5)  Dong;  L.;  Agarwal;  A.  K.;  Beebe;  D.  J.;  Jiang;  H.  R  Nature  2006; 
442,  551. 

(6)  Zhang;  Y.;  Li;  N.;  Delgado;  J.;  ZhoU;  N.;  Yoshida;  R;  Fraden;  S.; 
Epstein;  1.  R;  Xu;  B.  Soft  Matter  2012;  8,  7056.  Yuan;  P.  X.;  Kuksenok; 

O. ;  GrosS;  D.  E.;  BalazS;  A.  C.;  Moore;  J.  S.;  NuzzO;  R.  G.  Soft  Matter 
2013;  9;  1231. 

(7)  Zhang;  Y.;  ZhoU;  N.;  Akelk;  S.;  Kuang;  Y.;  Kim;  D.;  SchwartZ;  A; 
Bezpalko;  M.;  Foxman;  B.  M.;  Fraden;  S.;  Epstein;  1.  R;  et  al.  Angew. 
Chem.,  Int.  Ed.  2013;  52;  11494. 

(8)  Field;  R  J.;  NoyeS;  R.  M.;  KoroS;  E.  J.  Am.  Chem.  Soc.  1972;  94, 
8649.  Epstein;  1.  R;  Vanag;  V.  K.;  BalazS;  A.  C.;  Kuksenok;  O.;  Dayal; 

P. ;  Bhattacharya;  A.  Ace.  Chem.  Res.  2012,  45,  2160.  Bansagi;  T.;  Vanag; 
V.  K.;  Epstein;  1.  R  Science  2011;  331;  1309. 

(9)  Beebe;  D.  J.;  Moore;  J.  S.;  Bauer;  J.  M.;  Yu;  Q.;  LiU;  R  H.; 
DevadosS;  C.;  JO;  B.  H.  Nature  2000;  404,  588. 

(10)  Beebe;  D.  J.;  Mensing;  G.  A;  Walker;  G.  M.  Annu.  Rev.  Biomed. 
Eng.  2002,  4,  261.  Wu;  H.  K.;  Odom;  T.  W.;  ChiU;  D.  T.;  WhitesideS; 
G.  M.  /.  Am.  Chem.  Soc.  2003,  125,  554.  Zheng;  B.;  Roach;  L.  S.; 
IsmagiloV;  R  F.  /.  Am.  Chem.  Soc.  2003;  125,  11170.  Duffy;  D.  C.; 
McDonald;  J.  C.;  Schueller;  O.  J.  A;  WhitesideS;  G.  M.  Anal.  Chem. 
1998;  70,  4974. 

(11)  KarnS;  K.;  Vogan;  J.  M.;  Qin;  Q.;  Hickey  S.  F.;  Wilson;  S.  C.; 
Hammond;  M.  C.;  Herr;  A.  E.  /.  Am.  Chem.  Soc.  2013;  135;  3136. 
Puigmarti-LuiS;  J.;  Rubio-MartineZ;  M.;  Hartfelder;  U.;  ImaZ;  L; 
Maspoch;  D.;  Dittrich;  P.  S.  /.  Am.  Chem.  Soc.  2011,  133,  4216.  Price; 
G.  M.;  Chu;  K.  K.;  TrusloW;  J.  G.;  Tang-Schomer;  M.  D.;  Golden;  A. 
P.;  MertZ;  J.;  Tien;  J.  /.  Am.  Chem.  Soc.  2008,  130,  6664. 

(12)  Rossow;  T.;  Heyman;  J.  A;  Ehrlicher;  A.  J.;  Langhoff;  A;  WeitZ; 
D.  A;  Haag;  R;  Seiffert;  S.  /.  Am.  Chem.  Soc.  2012,  134,  4983. 

(13)  He;  M.;  Novak;  J.;  Julian;  B.  A;  Herr;  A.  E.  /.  Am.  Chem.  Soc. 

2011,  133,  19610.  He;  M.;  Herr;  A.  E.  /.  Am.  Chem.  Soc.  2010;  132; 
2512. 

(14)  Zhang;  H.;  Tumarkin;  E.;  Peerani;  R;  Nie;  Z.;  Sullan;  R.  M.  A; 
Walker;  G.  C.;  Kumacheva;  E.  /.  Am.  Chem.  Soc.  2006;  128,  12205. 

(15)  Wu;  H.  K.;  Huang;  B.;  Zare;  R  N.  /.  Am.  Chem.  Soc.  2006,  128, 
4194. 

(16)  Runyon;  M.  K.;  Johnson-Kerner;  B.  L.;  Ismagiloy  R  F.  Angew. 
Chem.  Int.  Ed.  2004,  43,  1531. 

(17)  Toiya;  M.;  Vanag;  V.  K.;  Epstein;  1.  R  Angew.  Chem.  Int.  Ed. 
2008;  47;  7753.  Delgado;  J.;  Li;  N.;  Leda;  M.;  Gonzalez-Ochoa;  H.  O.; 
Fraden;  S.;  Epstein;  1.  R.  Soft  Matter  2011,  7,  3155. 


7346 


dx.doi.org/10.1021/ja503665tU.  km.  Chem.  Soc.  2014,  136,  7341-7347 


Journal  of  the  American  Chemical  Society 


Article 


(18)  Ginri;  B.  T.;  Steinbock;  B.;  Kahveci;  M.;  Steinbock;  O.  /.  Phys. 
Chem.  A  2004,  108,  1325. 

(19)  Cheri;  1.  C.;  Kuksenok^  O.;  Yashin^  V.  V.;  Mosliri;  R.  M.;  BalazS; 
A.  C.;  Van  Vliet,  K.  J.  Soft  Matter  2011,  7,  3141. 

(20)  Xia;  Y.  N.;  Whitesides^  G.  M.  Angew.  Chem.  Int  Ed.  1998;  37, 
551. 

(21)  Park;  S.;  Kim;  D.;  Ko;  S.  Y.;  Park;  J.;  Akelk;  S.;  Xu;  B.;  Zhang;  Y.; 
Fraden;  S.  Lah  Chip  2014;  14,  1551. 

(22)  Chilkoti;  A.;  Dreher;  M.  R.;  Meyer;  D.  E.;  Rancher;  D.  Adv.  Drug 
Delivery  Rev.  2002,  S4,  613.  Sershen;  S.  R.;  Westcott;  S.  L.;  HalaS;  N.  J.; 
West;  J.  L.  /.  Biomed.  Mater.  Res.  2000;  SI,  293.  Stile;  R.  A;  Burghardt; 
W.  R.;  Healy;  K.  E.  Macromolecules  1999,  32,  7370. 

(23)  Delgado;  J.;  Zhang;  Y.;  Xu;  B.;  Epstein;  1.  R.  /.  Phys.  Chem.  A 

2011,  ns,  2208. 

(24)  Burstall;  F.  H.  /.  Chem.  Soc.  1936,  176. 

(25)  Cho;  S.;  Ha;  S.  H.  Struct.  Multidiscip.  Optim.  2009,  38,  53. 

(26)  Vanag;  V.  K.;  Epstein;  I.  R.  Phys.  Rev.  E  2011,  84.  Toiya;  M.; 
Gonzalez-Ochoa;  H.  O.;  Vanag;  V.  K.;  Fraden;  S.;  Epstein;  I.  R.  /.  Phys. 
Chem.  Lett.  2010,  1,  1241.  Taylor;  A.  F.;  Tinsley;  M.  R.;  Wang;  F.; 
Huang;  Z.  Y.;  Showalter;  K.  Science  2009;  323,  614. 

(27)  TompkinS;  N.;  Li;  N.;  Girabawe;  C.;  Heymann;  M.;  Ermentrout; 
G.  B.;  Epstein;  1.  R.;  Fraden;  S.  Proc.  Natl.  Acad.  Sci.  U.  S.  A.  2013;  111, 
4397. 

(28)  Eddington;  D.  T.;  LiU;  R.  H.;  Moore;  J.  S.;  Beebe;  D.  J.  Lab  Chip 
2001,  1,  96. 

(29)  Aihara;  R.;  Yoshikawa;  K.  /.  Phys.  Chem.  A  2001;  lOS,  8445. 


7347 


dx.doi.org/10.1021/ja503665tU.  Am.  Chem.  Soc.  2014,  136,  7341-7347 


